A simple model for semipermeable membrane: Donnan equilibrium 
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We study a model for macroions in an electrolyte solution confined by a semipermeable membrane. 
The membrane finite thickness is considered and both membrane surfaces are uniformly charged. 
The model explicitly includes electrostatic and size particles correlations. Our study is focused on 
the adsorption of macroions on the membrane surface and on the osmotic pressure. The theoretical 
prediction for the osmotic pressure shows a good agreement with experimental results. 
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I. INTRODUCTION 

Physics of two ionic solutions separated by a semiper- 
meable membrane is of wide interest in cell biology and 
colloids science P, 0, Q. A thermodynamical study 
of this problem was first carried out by F. G. Donnan 
0, considering the two fluids phases (here referred 
as a and j3) in the following way: (i) the a-phase con- 
tains two (small) ionic species, (ii) the /3-phase contains 
the same ionic species as the a-phase plus one macroion 
species. The two fluid phases are separated by a mem- 
brane which is permeable to the small ions and imperme- 
able to macroions, therefore, they interchange small ions 
whereas macroions are restricted to the /3-phase. The 
permeability condition is imposed by assuming (in both 
phases) a constant chemical potential of the permeating 
species. In this simple model, Donnan derived an expres- 
sion for the osmotic pressure (in terms of the ionic charge, 
concentration and excluded volume) which well describes 
systems close to ideality. Historically, this problem has 
been known as Donnan equilibrium. More recently some 
theories have been proposed to interpret osmotic pressure 
data 0, . These theories consider phenomenologically 
the macroion-macroion and ion-macroion many-body in- 
teractions, and provide a better fit for the osmotic pres- 
sure of macroions solutions than Donnan theory. 

The surface of a biological membrane has a net charge 
when it is in aqueous solution Q thus, at a fluid- 
membrane interface a broad variety of phenomena oc- 
cur. It is known that ionic solutions in the neighbor- 
hood of a charged surface produce an exponentially de- 
caying charge distribution, known as the electrical dou- 
ble layer. Low-concentrated solutions of monovalent 
ions are well described by the Gouy-Chapman theory 
(Poisson-Boltzmann equation) 8, 9]. However, multiva- 
lent ions display important deviations from this picture 
|lCl|| and more powerful theories from modern statisti- 
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cal mechanics (such as molecular simulation 0, 0, ^| , 
density functional 0, 0, 0| and integral equations 
EllllllliliJl) have been implemented for studying 
ions adsorption and interfacial phenomena. Macroions 
adsorption is a subject of current interest: In molecu- 
lar engineering, the macroions adsorption mechanisms 
are basic in self-assembling polyelectrolyte layers on a 
charged substrate [23| and novel colloids stabilization 
mechanisms j23l |. 

By means of integral equations, a previous study of 
Donnan equilibrium has been carried out by Zhou and 
Stell T hey used the method proposed by 

Henderson et al. which can be described as fol- 

lows: starting from a semipermeable spherical cavity 
|26l |. the planar membrane is obtained taking the limit 
of infinite cavity radius. Within this model, they ob- 
tained the charge distribution and mean electrostatic 
potential. However, due to the approximations used, 
they end up just with the integral version of the lin- 
ear Poisson-Boltzmann equation. A general shortcoming 
of the Poisson-Boltzmann equation is that ionic size ef- 
fects (short range correlations) are completely neglected, 
in consequence, the description of interfacial phenomena 
and computation of thermodynamical properties is lim- 
ited and valid only for low values of charge and concen- 
tration. 

From previous studies of two fluid phases separated 
by a permeable membrane, it is known that the adsorp- 
tion phenomena are strongly influenced by the mem- 
brane thickness [23, l28L l29j . On the other hand, short 
range correlations influence effective colloid-colloid inter- 
action [33, 1^ 113, 113, 13 , thermod yna mical properties 
35] and adsorption phenomena |3a . l37| in colloidal dis- 
persions. From these antecedents it is seen that there 
are several relevant aspects not considered in previous 
studies of Donnan equilibrium which deserve a proper 
consideration. In this study we consider explicitly the 
following effects: many-body (short and long range) cor- 
relations, the membrane thickness and the surface charge 
densities on each of the membrane faces. We use sim- 
ple model interactions and our study is carried out by 
means of integral equations. The theory gives the parti- 
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cles distribution in the neighborhood of the membrane, 
from which, the osmotic pressure is calculated. There 
are two points that we will address in this study: the ad- 
sorption of macroions at the membrane surface and the 
computation of the osmotic pressure for macroions solu- 
tions. Concerning the adsorption phenomena, we observe 
a broad variety of phenomena: charge reversal, charge in- 
version and macroions adsorption on a like-charged sur- 
face due to the fluid-fluid correlation. The computed 
osmotic pressure is compared with experimental results 
for a protein solution, obtaining an excellent agreement 
over a wide regime of concentrations. 

The paper is organized as follows: In section ^] we 
describe the integral equations method and the mem- 
brane and fluid models. In the same section, we derive 
the hypernetted chain/mean spherical (HNC/MS) inte- 
gral equations for the semipermeable membrane and the 
equations to compute the osmotic pressure. In section 
lllll a variety of results are discussed and finally in section 
II VI some conclusions are presented. 



II. THEORY 
A. Integral equations for inhomogeneous fluids 

The method that we use to derive integral equations 
for inhomogeneous fluids makes use of a simple fact: In a 
fluid, an external field can be considered as a particle in 
the fluid, i. e., as one more species infinitely dilute. This 
statement is valid in general, however, it is particularly 
useful in the statistical mechanics theory for inhomoge- 
neous fluids 0, described below. 

The multi-component Ornstein-Zernike equation for a 
fluid made up of n -f 1 species is 

n+l „ 

%(r2i) = Cy(r2i)+^p,„ hi„i{r23)cmj{ri3)dv3, 

m=l 

(1) 

where pm is the number density of species m, hij(r2i) = 
5ij(r2i) — 1 and Cy (r2i) are the total and direct corre- 
lation functions for two particles at r2 and ri of species 
i and j, respectively; with gij{r2i) the pair distribution 
and r2i — r2 — ri. Among the most known closures 
between hij{r2i) and Cij{r2i) used to solve Eq. we 
have|3^ 

Cy(r2i) = -/3Mii(r2i) + /i.y(r2i) - ln(/y(r2i), (2) 
Cjj(r2i) = /y(r2i)exp{-/3u,j(r2i)}.g.y(r2i), (3) 
Cij(r2i) = -(3uij{r2i) for r2i = |r2i | > . (4) 

Eqs. Q to Q are known as the hypernetted chain 
(HNC), the Percus-Yevick (PY) and the mean spherical 
(MS) approximations, respectively; My(r2i) is the direct 
interaction potential between two particles of species i 
and j, aij is their closest approach distance, fij (r2i) = 



exp{— (r2i)} — 1, and /? = l/ksT. Some more pos- 
sibilities to solve Eq. are originated by considering 
a closure for Cij(r2i) in the first term of Eq. and a 
different one for c„y (ris) in the second term of Eq. 
giving rise to hybrid closures. 

To derive integral equations for inhomogeneous fluids, 
we let an external field to be one of the fiuid species, 
say (n -|- l)-species (denoted as the 7-species), which is 
required to be infinitely dilute, i. e., 0. Therefore, 

the total correlation function between a 7-species particle 
and a j-species particle is given by 

n „ 

/J7j(r2i) = c^j(r2i) + ^ Pm hj„i{r23)c„ij{ri3)dv3 

m=l 

with j = 1, 71.(5) 

The total correlation functions for the remaining species 
satisfy a n-component Ornstein-Zernike equation as 
Eq. JU) (with no 7 species) from which c,„j(ri3) is ob- 
tained. In this scheme, the pair correlation functions, 
571 (r2i), is just the inhomogeneous one-particle distribu- 
tion function, gj(ri), for particles of species j under the 
influence of an external fleld. Thus, hjj{r2i) and Cjj{r2i) 
can be replaced with /ij(ri) = .gj(ri) — 1 and Cj{ri), re- 
spectively. Thus, the inhomogeneous local concentration 
for the j species is given by 

Pj(i"i) = Pj5j(ri), (6) 

By using the HNC closure (Eq. Q) for Cjj{r2i) in 
Eq. lO, we get 

gj(ri) = exp |-/3uj(ri) + ^ p„ ^ /im(r3)cmj(ri3) j 

(7) 

where the subindex 7 has been omitted for consistency 
with Eq. ®. In our approach, Cmj(ri3) in the inte- 
gral of Eq. (O, is approximated by the direct correlation 
function for a n-component homogeneous fluid. Thus, 
Cmj (ris) is obtained from Eq. |^ using one of the closures 
provided by Eqs. ©-Q- For the present derivation we 
will use Cmj {r 13) obtained with the MS closere (Eq. Q), 
therefore, we obtain the hypernetted chain/mean spher- 
ical (HNC/MS) integral equations for an inhomogeneous 
fluid. This equation has shown to be particularly success- 
ful in the case of inhomogeneous charged fluids when it is 
compared with molecular simulation data 0, 013^3 ■ 



B. The semipermeable membrane and fluid models 

The membrane is modelled as a planar hard wall of 
thickness d, and charge densities ai and (T2 on each sur- 
face and separates two fluid phases, referred as a and 
(3. In Fig. n (Ti and the a-phase are at the left hand 
side, whereas 172 and the /J-phase are at the right hand 
side. The fluid phased are made up in the following way: 
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I 1 particle 3, species m 




particle 1, species j 




particle 2, species 7 



being pf and pf the bulk concentrations of the i species 
in the a and (3 phases, respectively. The charge on the 
membrane is compensated by an excess of charge in the 
fluid (per unit area) , a' [3 Oil : 



cr^ = -(JT, 



(11) 



with (It = CTi + (72 and being cr" and the excess of 
charge in the a-phase and /3-phase, respectively, which 
are given by 



and 



being 



Pel{x)dx 



Pel{x)dx, 



(12) 



(13) 



FIG. 1: Schematic representation for a model macroions solu- 
tion confined by a semipermeable membrane. For the integra- 
tion of Eq. H18^ it is considered the axial symmetry around 
the ri vector where we fix a cylindrical coordinates system 



(i) The a-phase is a two component electrolyte, consid- 
ering the components to be hard spheres of diameter 
with a centered point charge qi = ZiC (being Zi and e the 
ionic valence and the proton's charge, respectively) and 
i — 1,2 standing for the species number. The solvent 
is considered as a uniform medium of dielectric constant 
£. (ii) The /3-phase is considered in the same way as the 
a-phase (same solvent dielectric constant and containing 
the same ionic species) plus one more species (species 3) 
of diameter 03 and charge . For simplicity, we have con- 
sidered that the membrane dielectric constant is equal to 
that of the solvent. In addition it is considered that 



a = fli = 02 < 03, 

2=2:1= -22- 



(8) 



Hence, we refer to the third species as the macroion 
species. Two ions of species m and j with relative posi- 
tion r, interact via the following potential 



Umjir) 



00 for r < amj, 

.......2 (9) 



er 



for r > ami, 



with m,j — 1, . . . , 3 and amj = ^ — -■ Far away from 

the membrane each phase is homogeneous and neutral, 
thus the neutrality condition is written as 



(10) 



Pel{x) 



m—1 



(14) 



the local charge density profile. In Eqs. H12|l . H13|) 
and (|14() it has been considered that the particles re- 
duced concentration profile {gj{ri)) depends only on 
the perpendicular position to the membrane, x, i. e., 
gji^i) — gj{x). Hence, the local concentration profile is 

Pniix) = Pmgm{x). 

According to the integral equations method outlined 
in section FlI Al the membrane is considered as a the fluid 
species labelled as species 7 (see Fig.^. The interaction 
potential between the membrane and a j-species par- 
ticle depends only on the particle position, x, referred 
to a coordinates system set in the middle of the mem- 
brane and measured perpendicularly. Thus, we write 
Mj(ri) = Uj{x) which is split as Uj{x) — m|'(x) -I- u*{x), 
being w^'(a;) the direct electrostatic potential and u*{x) 
the hard-core interaction. The former can be found from 
Gauss' law, resulting 



-Puf{x) = 



^ZjeParix- L) 



for a; > |, 



^Zje(3 [(Tt{-x -L)- (cTi - a2)d] for x<^, 

(15) 



where L is the location of a reference point. The hard- 
core interaction is given by 




u*{x) 



for i — 1, 2. For the impermeable species {j = 3) 

d + 03 

I 00 lor X < 

ul{x) = 



(16) 




(17) 
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as 



This potential imposes gsix) = for x < 

In Eq. O we use the expression of Cmj(r 13) for a prim- 
itive model bulk electrolyte, which has an analytical ex- 
pression, written as 



with 8(a;) the step function, defined as 

„ . . f for a; < 0, 
^^^> ^11 for a; > 0. 



(24) 



ins) 



for ri3 > a„j, 

er 

{ri-i) + cjjy (ria) for ria < a^j, 

(18) 

where ria = jriaj is the relative distance between two ions 
of species m and j . The particles short range correlations 
are considered through the c^j(''i3) and c^''j(r 13). The 
explicit form of these functions is given in appendix A. 
The integral in Eq. ((TJ can be expressed in cylindrical 
coordinates, and analytically calculated in the </> and r 
variables (see Fig. i. e., we consider a cylindrical 



The expressions for the kernels, Kmj{x, y) and Dmj{x, y), 
are given in appendix A. 

From the solution of Eq. 1)19(1 one obtains the reduced 
concentration profile, Pj{x) = pjgj{x). The bulk concen- 
trations of species j, at the a and /? phases (p" and p^) 
are given by 



and 



Pj = lim P]9j{x), 



p" = lim Pjgj{x), 



coordinates system where rf^ 



y — 2xy and 



X — > — 00 



(25) 



(26) 



dv3 = ddrdrdy. After a lengthy algebra, from Eq. ((TJ we 

get [Mil 

f 27r 

gj{x) — exp < — ZjeP [ai + <J2)\x\ — 2n A j{x) 

2 ^- d+ 

+ 27r^p, 



rn— 1 

3 



rn^l — 2 — 



""m 

{y)Gm]{x,y)dy 
{y)Gra3{x,y)dy 



respectively. At the /3 phase gj{x) — > 1 as 2: — > 00 then 
Pj ~ pj. On the other hand, in the a phase (for x < 
0), lim gj{x) ^ 1- It must be pointed out that the 

electrolyte bulk concentration at the a-phase (p" , for j = 
1, 2) satisfy the bulk electroneutrality condition, Eq. I|1U|1 . 
and are a result from the theory. 



(19) 



27rz. 



27r0. 



^ ^ ^771 Pn 



E 

m— 1 



9yn{y)[y+ \x-y\]dy 



hm{y)[y +\x~ y\]dy 



C. Computation of the osmotic pressure: contact 
theorem 

Let us consider a slice of fluid of width dx, area of its 
faces. A, parallel to the membrane and located at x. The 
force on the slice in the x direction, dF^, is given by 



The first and third integrals include hj (y) — gj {y) — 1 for 
particles in the a-phase whereas, the second and fourth 
integrals, for particles in the /3-phase. Notice the different 
summation limits due to the different phase composition. 
We have defined 



dF^{x) = E^{x)dQ + Sdp{x) 



(27) 



being Ex{x) the electric field in the x direction at x, dQ 
the total charge in the fluid slice and dp{x) the pressure 
difference at the two faces. Taking into account that 



Gmj{x,y) - 
Lmj{x,y) = 

K,nj{x,y) = 
and 

Ajix) 



LTnj{x,y) +Kmj{x,y), 



cZj{'>^i3)ri3dri 



c'^iiri3)ri3dr 



13, 



(20) 

D^j{x,m) 



(22) 



Ex{x) 



dij}{x) 
dx ' 



and using Poisson's equation, 

d'^i^ix) 47r 



dx'^ 



Pel {x), 



(28) 



(29) 



we write 



P3 
2 

771—1 



G3j{x,y)dy 



dQ = Ap,i{x)dx^-^(^^P^]dx. (30) 



Thus, we rewrite Eq. (|27|l as 



47r V dx^ 



Grnjix,y)dy 



+ z. 



/3e' 



2:3^3 — / {y 



dF,{x) = — ( ] l-^ ] dx + Adp{x), (31) 



y\]dy 



At: \ dx 
or equivalently, 



ZiC- 



(3d 



(ai -a2)e(x-f d/2), 



(23) 



dF^ix) = — TT 



As d f dif}{x) 



n dx \ dx 



dx + Adp{x). (32) 
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Considering that in equilibrium dF^ {x) = and integrat- 
ing Eq. (|32|1 in the interval [d/2,oo) with the boundary 
condition 



dtp{x) 



lim 



0, 



we obtain 



£ / dip{x) 
8^ V dx 



(33) 



(34) 



where 11'' = limj^^aoPix) is the bulk fluid pressure and 
the expression for the pressure on the membrane right 
surface, = p{0), is given by 



Po 



^kBTp^{0) = kBTY, 



Pi9i 



(35) 



where p^(0) = X^Li PiC^)- Eq. lESJ is an exact rela- 
tionship which can be obtained by considering the force 
on the fluid (at the contact plane) exerted by the hard 
wall . From basic electrostatics we have 



e / dtp{x) 
An 



dx 



Pel {x)dx, 



=d/2 •'d/2 

thus, using Eqs. H34|) . H35|) and (|36|l we can write 



(36) 



d/2 



Pel {x)dx 



ksTy^Pigi 

1=1 



(37) 

where the first term can be identified as the Maxwell 
stress tensor. A similar expression is obtained for the 
bulk pressure in the a phase 



e 



-d/2 



Pel{x)dx 



(38) 



The osmotic pressure, 11, is defined as 

n = n'' - n". (39) 

From of Eqs. ^ and Eq. ^ becomes, 



Pi9i 

1=1 



a. 



(40) 



Due to the fluid-fluid correlation across a thin mem- 
brane, the induced charge densities((T" and a^) and 
gi (±^i2ij depend on the fluid conditions {pt, Zi, at, with 
i = 1, 3) and on the membrane parameters (cti, (72 and 
d) [H 121 113. The computed value of H (using a", 
and gi {±^^) from HNC/MS), however, does not de- 
pend on the membrane parameters. We have numerically 




FIG. 2: Reduced concentration profiles (RCPs) for a 
macroions solution {pM = O.OIM, zm ~ —10) in a monova- 
lent electrolyte (p+ = I.IM and p_ = l.OM), with 03 = 3.8a, 
o"i = o"2 = 0.272 C/m^and d = a. The continuous, dashed 
and dotted lines represent the RCPs for the macroions, anions 
and cations, respectively. 



corroborated this fact by computing 11 for several values 
of ai, (72 and d. This is physically appealing since the 
pressure can only depend on the bulk fluid conditions at 
both sides of the membrane. However, we had to do very 
precise calculations of gi{x), particularly in the neighbor- 
hood of a; = ± , to prove the above statement. 

Eq. (|4()|l is an exact theorem to compute the osmotic 
pressure, 11, in terms of microscopic quantities. A sim- 
ilar expression for the osmotic pressure was derived by 
Zhou and Stell [l^. However, the differences between 
the current derivation and that of those authors are due 
to the ions-membrane short range interactions. If we use, 
in the Zhou and Stell theory, the hard wall interaction 
between the permeable ions and the membrane (provided 
by Eq. ^fTE\i ) we recover Eq. l|in|) . 



III. RESULTS AND DISCUSSION 

Several physical effects determine particles adsorption 
on the charged membrane. One of the most relevant is 
the membrane-particle direct interaction energy, which, 
at the surface is given by 



(41) 



being L the location of a reference point. The more 
negative the value of C/^ < 0, particles adsorption is 
energetically more favorable. Many body correlations 
play also an important role in the adsorption phenomena 
and are responsible for the surface-particle forces of non- 
electrostatic origin. Although it is not possible to sharply 
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x[a/2] x[a/2] 



FIG. 3: Same as in Fig.|5|but with 03 = 7a. The hnes meaning 
is the same as in Fig. |21 



distinguish the origin of correlations, we assume that the 
particles volume fraction {r]T = f X^i Pi^-i) quantifies the 
contribution of short range correlations. To quantify the 
effect of the coulombic interaction (long range correla- 
tions), we define the parameter = I3qiqj/eaij and we 
simply write = = /Sqf/eai, when i = j. In the dis- 
cussion, the effects produced by varying the fluid param- 
eters (ai, Zi and pi, with i = 1,...,3) will be associated 
with the increment (or decrement) of the contributions 
arising from short and long range correlations {rjT and 

For our discussion it is useful remark the following 
concepts: When the amount of adsorbed charge exceeds 
what is required to screen the surface charge, it is said 
to occur a surface charge reversal (CR). In consequence, 
at certain distance from the surface the electric field is 
inverted. Next to the CR layer, a second layer of ions 
(with same charge sign as that of the surface) is formed, 
producin g a charge inversion (CI) of the electrical double 
layer [iM l49j . Such a denomination (charge inversion) is 
originated from the fact that ions invert their role in the 
diffuse layer. In the past, we have shown that charge 
reversal and charge inversion are many body effects and 
are induced by a compromise between short range cor- 
relations (^t) and electrostatic long range correlations 

The HNC/MS equations for a semipermeable mem- 
brane, Eqs. (|19|l , are numerically solved using a finite 
element technic |5l|,|53]. From the solution of HNC/MS 
equations the reduced concentration profiles (RCPs), 
gj{x), are obtained. In the discussion, we adopted the 
following notation: zi = z^, Z2 — z~ and 23 — zm for 
the number valence of cations, anions and macroions, 
respectively; idem for pi and a^. In all our calcula- 
tions we have used fixed values of T = 298K, e = 78.5 
and a = 4.25 A. The effects of salt valence, macroions 



FIG. 4: Reduced concentration profiles (RCPs) for a 
macroions solution (pA/ = O.OIM, zm ~ —10) in a monova- 
lent electrolyte (p+ = I.IM and p_ — l.OM), with 03 — 3.8a, 
ai = a2 ~ —0.272 C/m'^and d — a. The lines meaning is the 
same as in Fig. |5| 



size, and membrane surface charge density on macroions 
adsorption are analyzed. Thus, we consider macroions 
(pa/ — O.OIM and zm — 10) in a (a) monovalent (p+ = 
I.IM and p_ = l.OM) and (b) divalent (p+ = 0.55M and 
p_ = 0.5M) electrolyte solutions. For each case, we con- 
sidered two macroion sizes (qm = 3.8a and gm = 7a) 
and three values for the membrane charge densities: i) 
(Ti = C72 ^ 0.272C/m2, ii) cri (72 = -0.272C/m2 and 
iii) cTi = 0.68C/m2, a2 = -O.lSGC/m^. Finally we show 
calculations for the osmotic pressure (as a function of the 
macroion concentration) compared with experimental re- 
sults. 



A. Macroions in a monovalent electrolyte {z = 1) 

In this subsection we discuss the case of macroions 
(zm = —10) in a monovalent electrolyte solution (p+ = 
I.IM, p_ = l.OM, z+ = -z = 1). 



1. Positively charged membrane 

In Fig. 121 we present the RCPs when the membrane 
is positive and symmetrically charged {ai =02 = 
0.272C/m~), the membrane thickness is c? = a and 
o-M = 3.8a. As a result, we obtained the asymp- 
totic values of the distribution functions in the a-phase, 
g-{-oo) = 1.0758 and 5+(-oo) = 0.9779, such that 
Eq. (|10ll is satisfied. In both phases we observe that 
the negative ions are adsorbed and the positive ions ex- 
pelled, as it is expected. In the /3-phase we observe that 
the adsorption of macroions is more favorable than the 
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adsorption of negative small ions. This is understood in 
terms of the electrostatic energy of one particle of species 
i at the membrane surface, Ui: since \zm\ > from 
Eq. H41|l it is easy to see that Um < U- < 0, which favors 
macroions adsorption. As we pointed out above, many 
body correlations also influence adsorption. In this case 
macroion-macroion long range correlations are predomi- 
nantly more important than the ion-ion correlations since 
S^jxi Ki 7^_. Short range correlations are also important 
and play an important role in the macroions adsorption: 
an increment of macroions concentration (keeping con- 
stant gm and zm) implies an increment of rjT and pro- 
duces an increment of the macroions adsorption. How- 
ever an increment of rjT not always is followed by an 
increase of the adsorption as it will be noticed in the 
discussion of Fig. O In the macroions RCP we find a 
first macroions layer next to the membrane surface, af- 
ter this layer there is a region where the macroions are 
completely expelled and then a small second peak. The 
negative small ions are also adsorbed to the membrane's 
surface but their concentration is smaller than for the 
macroions. The RCP for positive ions show that these 
are completely expelled from the membrane right surface. 
However, a peak in the RCP is found at a; « 4.3a which 
is the ion-surface distance when there is a macroion in 
between. This peak implies an effective surface-cation 
attraction due to a field inversion caused by the surface 
CR. 

Fig. 13 shows the results obtained for the same sys- 
tem as in Fig. |21 except that the macroions are larger, 
Um — 7a. At the a-phase the distribution function have 
the same qualitative behavior as in Fig. 12 however, in 
this case g-{—oo) — 1.1938 and g+{—oo) = 1.0853, i. 
e., the amount of salt in the a-phase increases by in- 
creasing the macroions size. Respect to Fig. |2 a decre- 
ment of macroions adsorption is observed in Fig. O The 
increment of the macroions diameter implies a decre- 
ment of the macroions-macroions long range correlations 
(in this case ~ 2i^_) and an increment of rjT, and 
hence, an increment of the contributions arising from 
correlations of short range nature. Such an increment, 
however, does not increase (but decrease) the macroions 
adsorption. This is understood in terms of Eq. I|41() 
[UMia-M — 3.8a) < Um{o,m — 7a)] which implies that 
the adsorption is energetically less favorable for gm = 7 
than for smaller macroions with the same charge. 



2. Negatively charged membrane 

Fig.2]shows the RCPs when the membrane is negative 
and symmetrically charged (cri — a2 — — 0.272C/m^) 
for Um — 3.8a. The membrane thickness is d = a. 
The asymptotic value of the distributions function does 
not depend on the membrane's charge and thickness, 
hence, the gi{—oo) have the same value as in Fig. |21 
The RCP for small ions behave in a normal way in the 
sense that positive ions are attracted to the to the mem- 
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FIG. 5: Reduced concentration profiles (RCPs) for a 
macroions solution (pM = O.OIM, zm = —10) in a monova- 
lent electrolyte (p+ = I.IM and p- = l.OM), with as — 3.8a, 
(Ti = 0.68 C/m^, (72 = -0.136 C/m^and d = a. The lines 
meaning is the same as in Fig. |5| 



brane surface, whereas the negative ions are expelled 
from it. In the /3-phase, it is observed that macroions 
are expelled from the membrane surface for x < 6(a/2). 
After this zero concentration region, an small peak in 
the macroions RCP is observed at a; w 9 (a/2) indicat- 
ing an effective macroion-membrane attractive force and 
an slight surface CR produced by the small cations. At 
these electrolyte conditions (1:1 and p+ — IM with 
no macroions), cations display a monotonically decaying 
distribution profile hence, the oscillatory behavior 
of the RCPs for the small ions (in the /3-phase) is a con- 
sequence of the presence of macroions. By considering 
the macroions size qm = 7.0 (not shown), the qualita- 
tive behavior of the RCPs is similar to that of Fig. 01 
with the following differences: (i) the RCPs oscillations 
of the small ions at the /3-phase are of longer range and 
(ii) the macroions RCP maximum is higher. This fact 
points out the relevance of the effect of the particles size 
and concentration (particles volume fraction) in the effec- 
tive attraction between like charged particles in solution 

US mm El 13. 

From the analysis presented from Fig. |2] to 01 we see 
that macroions adsorption increases by increasing rjT, Cm 
and —Um- However, we point out the following findings: 
when macroions and the surface are oppositely charged, 
long range electrostatic correlations dominate over short 
range correlations. Hence, adsorption is enhanced by in- 
creasing Cm and/or —Um- On the other hand, when 
macroions and the surface are like charged the mecha- 
nism for macroions adsorption (mediated by the small 
ions) is mainly driven by short range correlations, thus, 
adsorption increases by increasing rjx even though Cm 
decreases. 
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FIG. 6: Induced charge densities as a function of the mem- 
brane thickness for a macroions solution (pM = O.OIM, 
zm = —10) in a monovalent electrolyte (p+ = I.IM and 
p- = l.OM) with (71 = 0.68 C/m^and (72 = -0.136 C/m^ 
The continuous line represents calculations for au = 3.8a, 
whereas the dashed line for om = 7.0o 



3. Unsymmetrically charged membrane 

In Fig. |5l it is shown the RCPs at the two membrane 
sides for a macroions diameter = 3.8a and d = a. The 
membrane is unsymmetrically charged with cti = 0.68 
C/m^and (72 = -0.136 C/m^. At the right hand side 
surface macroions and small anions (negatively charged) 
are adsorbed on the membrane, in spite of (72 < 0. The 
adsorption of negatively charged particles on a negatively 
charged surface is due to the correlation between the two 
fluids. This is understood by considering the following 
two facts: (i) the membrane has a positive net charge 
((7t > 0) and (ii) the charge on the left hand side surface 
((7i ) is not completely screened by the excess of charge in 
its corresponding fluid phase (cr"), therefore, the electric 
field produced by (7i + a" overcomes the field produced 
by a2, inducing macroions and anions adsorption. From 
here, we observed that cr" and depend on d as it is 
discussed below. 

In Fig. IHI it is shown the excess of charge densities 
(7q and ap, as a function of the membrane's thickness 
d. The dependence on the wall thickness of the induced 
charge densities cr^ and (7/3 , is a manifestation of the cor- 
relation between the fluids. The correlation between the 
two fluids is due to the electrostatic interaction among 
the particles at both phases but, more importantly, to 
the fact that they are at constant chemical potential. 
For a sufficiently large membrane's thickness the induced 
charge density in each fluid phase screens its correspond- 
ing membrane surface, i. e., (7" — > — (7i and —(72 as 
d ^ oo. At d = 100a each fluid has screened its respec- 
tive surface charge density. Here we show a comparison 
between the results obtained for ajvf = 3.8a and um — 7a. 
In both cases the results are qualitatively similar. 
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FIG. 7: Reduced concentration profiles (RCPs) for a 
macroions solution {pM = O.OIM, zm = —10) in a divalent 
electrolyte (p+ = 0.55M and p_ = 0.5M), with as = 3.8a, 
(7i = (72 = 0.272 C/m^and d = a. The continuous, dashed 
and dotted lines represent the RCPs for the macroions, anions 
and cations, respectively. 



B. Macroions in a divalent electrolyte (z = 2) 

We now discuss the case of macroions {zm = —10) in 
a divalent electrolyte solution (p+ = 0.55M, p_ = 0.5M, 

z+ = -z = 2). 



1. Positively charged membrane 

In Fig.Qwe show the RCPs for (7i = (72 = 0.272C/m2, 
d — a and ajv/ — 3.8a. At the (3!-phase we observe os- 
cillations of the RCPs, which is a typical behavior of a 
divalent electrolyte. Although we observe a strong ad- 
sorption of macroions in the /3-phase, the amount of ad- 
sorbed small negative ions is even larger (the concentra- 
tions of macroions and small negative ions at the inter- 
face are pM (^^^) « 7.2M and p_ {^) « 16M, re- 



spectively). Energetically, macroions adsorption should 
be more favorable, however, macroions adsorption is in- 
hibited because divalent positive ions more efficiently 
screen macroion-membrane and macroion-macroion in- 
teractions. We infer this from the RCP for cations, which 
displays two small peaks at x « 2a and x ~ 4.3a. The 
first peak corresponds to a positive ions layer contiguous 
to the negative ions adsorbed on the wall. The position 
of the second peak corresponds to a positive ions layer 
next to the macroions layer. This structure indicates that 
positive ions surround macroions due to their strong elec- 

2 

trostatic interaction, in this case £,m+ ~ sksTa ^ 
— 14, whereas for macroions in a monovalent electrolyte 
S.M+ ~ —7. By comparison of Figs. Inland [7| we see that 
macroions (next to an oppositely charged surface) are 
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FIG. 8: Same as in Fig.[7|but with 03 = 7a. The hnes meaning 
is the same as in Fig.Q 

better adsorbed when they are in a monovalent solution 
rather than in a multivalent solution: for this partic- 
ular case of macroions {pM = O.OIM, zm — —10 and 
aM = 3.8a), pm{^^-§^) = 20M when macroions are in 
a monovalent electrolyte, whereas pm ( ) « 7.2M 
when macroions are in a divalent electrolyte. 

In Fig.|Slwe show the RCPs for the same conditions as 
in Fig.|3but aM = 7a. Although macroions adsorb in the 
/3-phase, they do not influence significantly on the local 
concentration of small ions. Hence, we see that the RCPs 
for small ions are quantitatively similar in both phases: 
the concentrations of counterions at the membrane sur- 
faces are p (±^) « 22. 5M, in addition, the RCPs max- 
ima are located symmetrically around x ~ ±1.7a. The 
adsorption of macroions decrease (respect to Fig. Cj) due 
to the efficient screening of the membrane charge by the 
small negative ions and because the adsorption of larger 
ions (keeping zm constant) is energetically less favorable, 
as it was pointed out in the discussion of Fig. |21 We see 
that the layer of cations around macroions (seen in Fig.[7I) 
desapears due to the decrement of their coulombic inter- 
action, in this case Cm+ ~ —7.6. 

2. Negatively charged membrane 

In Fig. ini we show the RCPs for a negatively charged 
membrane (cti =172 = — 0.272C/m^) with d = a and 
o-M — 3.8a. The quantitative behavior of the small ions 
RCPs is similar at both phases, as much as the location of 
the maxima located symmetrically at x « ±1.8a. At the 
/3-phase, the macroions RCP displays a peak at x = 3.2a 
indicating the formation of a macroions layer. The ad- 
sorption of macroions at this layer is enhanced (respect to 
Fig.^ because of the surface charge rewersaZ produced by 
the divalent cations, i. e., macroions see the membrane 



FIG. 9: Reduced concentration profiles (RCPs) for a 
macroions solution (pM = O.OIM, zm = —10) in a divalent 
electrolyte (p+ = 0.55M and p_ — 0.5M), with 03 = 3.8a, 
(Ti = (72 = —0.272 C/m^and d — a. The lines meaning is the 
same as in Fig.[7| 

surface with an effective positive charge. By increasing 
macroions size (keeping zm constant), macroions attrac- 
tion to an oppositely charged surface is energetically less 
favorable, thus, we observe that macroions adsorption 
decreases by increasing macroions size, in opposition to 
the behavior of macroions in a monovalent electrolyte. 

It seems to be a general feature that macroions next to 
an oppositely charged surface are better adsorbed when 
they are in a monovalent solution than in a multiva- 
lent solution. It is important to point out the formation 
of a cations layer around the macroions when |^a/+ | is 
high (seen for divalent electrolyte in Fig. [7| ) . On the 
other hand, when macroions and the surface are like 
charged, a divalent electrolyte solution mediates an ef- 
fective membrane-macroions attraction, this attraction 
is favored by a higher value of rather than for a high 
value of rjM, as in the monovalent case. 



3. Unsymmetrically charged membrane 

In Fig. ^1 we show the RCPs at the two phase for 
(Ti = 0.68C/m2, (72 = -0.136C/m2, ag = 3.8a and 
d = a. The correlation between the two fluids is man- 
ifested by the attraction of negatively charged particles 
towards the negatively charged surface at the oi-phase. 
However, the adsorption of macroions is quite less effi- 
cient than in the monovalent electrolyte case of Fig. |31 
In the monovalent case the contact value of the local 
concentration is pm{^^^-§^) ~ 2.7M whereas in this case 
it is pm(^=^) ~ 0.12M. This is due mainly to the 
more efficient field screening by the divalent electrolyte 
at the left hand side surface. In the monovalent elec- 
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FIG. 10: Reduced concentration profiles (RCPs) for a 
macroions solution {pM = O.OIM, zm = —10) in a divalent 
electrolyte (p+ = 0.55M and p_ = 0.5M), with as = 3.8a, 
ai = 0.68 C/m^, a2 = -0.136 C/m^and d = a. The lines 
meaning is the same as in Fig. |7| 



FIG. 11: Reduced osmotic pressure as a function of the 
macroions concentration, pM, for zm = — 10 and two macroion 
sizes, aM = 5a and aM = 7a. The solid line represents the 
osmotic pressure for macroions in a monovalent electrolyte 
(p_ = l.OM) whereas the dashed line for macroions in a di- 
valent electrolyte (p_ = 0.5M). 



trolyte case the induced charge density at the a-phase 
is (J a = — 0.43C/ni^when d — a whereas in the divalent 
electrolyte case CTq, — — 0.50C/m^. 

C. The osmotic pressure 

Although the adsorption of macroions is strongly in- 
fluenced by the membrane surface charge and thickness, 
we observe however, that the osmotic pressure does not 
depends on these membrane properties as it has been 
pointed out in subsection III CI In Fig. ^2 we show the 
osmotic pressure (obtained from HNC/MS theory) as 
a function of the macroions concentration pM, for two 
macroion sizes. This plot shows the osmotic pressure for 
macroions in a monovalent electrolyte and in a divalent 
electrolyte. The osmotic pressure increases by increasing 
the particles excluded volume (either pu or om), on the 
other hand, it is not observed a qualitative difference be- 
tween the curves for the osmotic pressure of macroions 
in a monovalent and divalent electrolytes. 

Computation of the osmotic pressure of proteins so- 
lutions is an issue addressed by some authors 0, 113 ■ 
Results of particular interest are those for albumin so- 
lutions, for which theoretical calculations and measure- 
ments of the osmotic pressure have been reported. In 
Fig. 1121 we show the osmotic pressure predictions of 
HNC /MS theory (as a function of the protein concentra- 
tion) and experimental results. The experimental data 
correspond to a solution of albumin in a 0.15M NaCl 
aqueous solution. In accordance with titration measure- 
ments, the albumin has a net charge of Q = —9e and 
Q = — 20e for pHf* 5.4 and pHw 7.4, respectively, in 
our calculations we have used um = 62 A which corre- 



sponds to the experimental protein diameter. It is re- 
markable the excellent agreement between theory end 
experiment as well as the fact that no adjustable pa- 
rameters have been used. The prediction of HNC/MS 
fits well the experimental data even for protein con- 
centration as high as 30% the protein volume in solu- 
tion (which is estimated assuming the albumin molecu- 
lar weight Wai = 69Kg/mol). For higher protein concen- 
tration, HNC/MS shows discrepancies with experimen- 
tal measurements which may be associated with the fol- 
lowing facts: (i) the albumin molecule is not spherically 
symmetric as in the model, therefore, the protein geom- 
etry becomes relevant when the protein volume fraction 
is high, (ii) integral equations are approximated theories, 
meaning that they do not take into account all the par- 
ticle correlations. 



IV. CONCLUSIONS 

We studied a model charged membrane separating two 
fluid phases (a and /?). The /?-phase phase contains 
macroions in an electrolyte solution and the a-phase is 
a simple electrolyte solution. The system is modeled in 
such a way that the small ions at both phases are at 
the same chemical potential, thus, the membrane is con- 
sidered to be semipermeable. It is important to point 
out that we cionsidered explicitly the effect of parti- 
cles size (short range correlations) and electrostatic long 
range correlations. Here we have applied the hypernetted 
chain/mean spherical integral equations to the semiper- 
meable membrane model. Also we have derived the equa- 
tion for the osmotic pressure throug a simple forces bal- 
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FIG. 12: Bovine serum albumin reduced osmotic pressure as 
a function of albumin weight concentration, pMWai, at 25°C 
and in 0.15M NaCl at pH=7.4(Q = -20e) and pH=5.4(--9e) 
from ref. 0- The curves are obtained from Eq. 11211 using 
aM = 62A, a = 4.25A and 2 = 1. 



By solving the HNC/MS integral equations we ob- 
tained the particles concentration profiles which allowed 
us to study the adsorption of macroions on the mem- 
brane. We analyzed the influence of several factors in 
the macroions adsorption: membrane surface charge, 
membrane thickness, the effect of salt and macroions 
size. From this study we emphasize the following re- 
sults: 1) When the membrane and macroions are oppo- 
sitely charged the adsorption is energetically favorable, 
however, if multivalent ions are present the membrane- 
macroions interaction is screened which is unfavorable 
for macroions adsorption. On the other hand, the larger 
the macroions (keeping the same charge and concentra- 
tion) the adsorption is energetically less favorable. 2) 
The attraction of macroions towards a like-charged sur- 
face seems to be energetically unfavorable. Nevertheless, 
in our model, we find that such an attraction is feasible 
and it is due to short range correlations, which are prop- 
erly considered in our theory. These previous results are 
not a consequence of the permeability condition but are 
general. 3) As a consequence of the permeability condi- 
tion and constant chemical potential, we have that for 
an unsymmetrically charged membrane the fluids corre- 
lation may produce adsorption of charged macroions on 
a like charged surface. In addition, the permeability con- 
dition implies a non trivial relation between the induced 
charge densities (cr" and a^) and the membrane thick- 
ness. The theory predictions are robust as it is shown by 
the excellent agreement between theory and experiment 
where we have not used adjustable parameters. The re- 
sults of this work could be techn olog ically relevant for 
the design of selective membranes |54j . 
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APPENDIX A: APPENDIX 



1. The Mean Spherical Closure 



The primitive model is the simplest model for an elec- 
trolyte that includes many relevant aspects of real solu- 
tions. In the general case the primitive model is consti- 
tuted by n-species of particles, with the mixture is em- 
bedded in a uniform medium of dielectric constant e at 
temperature T . Each species is defined by the particles 
point charge at the center, qi = z^e (where e stands for 
the proton's charge and Zi for the ionic valence) , the ionic 
diameter, a^, and number concentration, pi. The fluid is 
constrained to the following condition 



ZiPi = 0. 



(Al) 



i=l 



The expressions for the direct correlation functions, 
Cy(ri3), for a feit/fc electrolyte (required in n Eq. ^) were 
obtained by Blum [s^ and Hiroike [s^, through the MS 
closure, and are written as 



c.,(ri3) = ^d.,(ri3) + 4^(ri3) - /3^, (A2) 

^ ^* 13 



with cf[(ri3) 



dijirii), (3 = l/ksT and 



+ forO<ri3<A.„ 



ri3 



- 5g) + 5gV,3 + 5gV?3, for A., < n3 < a, 
0, for ri3 > aij] 



with Xi 



and Oi 



fli + a," 



(A3) 

The constants 



12 



in Eq. (jA3|l are given by 



u(2) 



bf^ = {xi - Xj) (n, - rij) 



= (ni + Txi), 



[a^ -aj)\ [s, - Sjj 



(«» ~ «i) 
16 



[(jij + Vxi + Tij + Yxj)^ — AniUj] 



(jA7p are given by 

Ai = (l-77T)"^{l + ryT + ?7l + ^aVT[l + 2?7T] 

- 1^3(03 - a)^{a{l + 1J3) + a3[l + 2(771 + 772)]}} 

3 

+ ^i^-VryHpTi^ + VT + VT) (A8) 

TT TT ^ 1 

- -^^3(^1 +P2)(a3 - af[{a + a3)+aa3-^piaf] \ , 



5(4) 

b: 



(5) 



— s, + ^Sj+ n.Uj - - s- + , 
ai Gj Z ■' 

J * 



a = -71013313(013) Piaigiijai), 



1=1 



12 



(A9) 



(AlO) 



Bi = B2 = -TT [(pi +p2)a^gn(a) +P3a35?3(ai3)] , (All) 
with 



where a^^ are defined as Xi = Zi + riitti and F is obtained Snl*^) — .922(a) 
from the solution of the following algebraic equation 

ffi3(ai3) 



+ -^"^sal (o - 03) \ (1 - ?7t)~^ , 



1 = > Pi(Zj + niOi) . 



(A4) 



The rii are obtained from the solution of the following 
set of algebraic equations 



[asgiiia] + Qg33(a3)] 
2oi3 



(A12) 



The expressions for ^3, ,63 and 333 (03) are obtained by 
interchanging 771+772, pi + P2 and oi with 7^3, and 03, 
respectively, in the expressions for Ai Bi, gii(a). 



2. The kernels expressions 



{zi + 7^iOi) r = 7ii + coj^ (zj + 7ijaj), (A5) 



where c= -[1 - - E Pi4]~^- 
z bj=i 

Considering that = 01=02, c^j''(7'i3) is just the 
direct correlation function for a hard spheres binary mix- 
ture in the PY ap prox imation. For particles of the same 
size it is given by |53| 



-A^ - BiVis - Srf^ for ri3 < 



for ri3 > ai 



(A6) 



For particles of different size we have 

'Ai for s < Ai3, 



Carrying out the integrations indicated in Eqs. H21|l 
and (I22Jl, using Eqs. (|A3|) . (|A6|) and HA7|) . the expressions 
for Kij{x,y) and Dij{x,y) are 

D,j{x,y) = 

f bg^fco + Ji + b^f Afi ~ bgWz + bfMs + b[fM5, 
for < \x - y\ < Xij, 



ib^ + z,z,)j,-i^;>j, + b^>j, + b)^>j, 

for A.y < \x -y\ < aij, 
, 0, for a^j < \x -y\ 



(4) 



r,(5) 



(A13) 



^Ki,{x,y) 



A; J2 + J3 + 5 J5, for O;., > - ?/|, 



0, 



for Oij < |x - 7/1, 

(A14) 



C?|(n3) 



-Al 







for Ai3 < ri3 < oi3, 
for ri3 > 013. 

(A7) 



with X = ri3 — A13. The constants used in Eqs. ljA6|l and 



-Ki_z{x,y) = 

Ai J2 + ao^/a + (5Ai3a'* + (5a^/5, for \x - y\ < A13, 
Ai J2 + "^^3 + 4(5Ai3P4 + SP,^, for Ai3 <\x-y\< 013, 
0, for 013 < \x-y\, 

(A15) 

where we use the following definitions: 

J„ = {a^ -\x^ 77I") /71, (A16) 

P„ = (a"-(|a;-y|-Ayr)/«, (A17) 
Mn = {a^^)^)/n, (A18) 



13 



and 



fco = (A?- -{x- yf)/2. 



(A19) 
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